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Abstract 

We consider an integrable conformally invariant two dimensional model associated 
to the affine Kac- Moody algebra s[3(C). It possesses four scalar fields and six Dirac 
spinors. The theory does not possesses a local Lagrangian since the spinor equations 
of motion present interaction terms which are bilinear in the spinors. There exists a 
submodel presenting an equivalence between a U{1) vector current and a topological 
current, which leads to a confinement of the spinors inside the solitons. We calculate the 
one-soliton and two-soliton solutions using a procedure which is a hybrid of the dressing 
and Hirota methods. The soliton masses and time delays due to the soliton interactions 
are also calculated. We give a computer program to calculate the soliton solutions. 



1 Introduction 



The non-perturbative aspects of Lorentz invariant field theories are related in one or another 
way to classical soliton solutions. The best known example is provided by the strong coupling 
sector of non-abelian gauge theories. It is believed that the fundamental particles, namely 
gluons and quarks, get confined and form bound states which should correspond to the known 
spectrum of hadrons. The most popular mechanism to explain such confinement is based on 
a condensation of magnetic monopoles producing a vacuum which is a magnetic supercon- 
ductor. The quarks and gluons get then confined by the dual Meissner effect. The monopoles 
responsible for that mechanism are not excitations of the fundamental fields appearing in 
the Lagrangian. They are classical soliton solutions carrying a topological charge. There are 
many ideas based on Montonen-Olive duality conjecture [0] where the monopoles are inter- 
preted as the fundamental particles of a dual Lagrangian describing the strong coupling sector 
of gauge theories. 

The development of exact methods for studying solitons and non-perturbative aspects 
of field theories is important for many reasons. Besides its relevance for gauge theories, it 
finds applications in many areas of condensed matter physics, non-linear phenomena and 
fluid dynamics. In many cases, the symmetries involved allow the reduction of the number of 
dimensions, and integrable non-linear models in low dimensions play an important role. The 
soliton theory in two dimensions has reached a high degree of development. It is known how to 
construct practically all the soliton solutions of two dimensional models using some very basic 
methods. Such theories have a zero curvature representation with the potentials taking values 
on an affine Kac-Moody algebra. In addition, it should present some "vacuum" solutions 
such that the zero curvature potentials belong to an abelian subalgebra when evaluated on 
them. The soliton solutions are then constructed via the dressing method @] where the 
dressing transformation are generated by the operators which diagonalize the adjoint action 
of the "oscillators" of that abelian subalgebra 0. Some of the structures of two-dimensional 
integrable models can in fact be generalized to theories in higher dimensions |Q. 

In this paper we consider a two-dimensional field theory involving four scalar fields and six 
Dirac spinors. The interactions among the fields present some peculiar properties. There are 
the usual coupling of bilinears in the spinors to exponentials of the scalars. However, some of 
the equations of motion for the spinors present terms which are bilinear in the spinors them- 
selves. That fact makes it difficult to believe that one can find a local Lagrangian for the 
theory. In spite of that fact, the model presents a lot of symmetries. It is conformally invari- 
ant, possesses local gauge symetries as well as vector and axial conserved currents bilinear in 
the spinors. In addition, it presents three species of one-soliton solutions carrying non-trivial 
topological charges. One of the most striking properties of the model is that it presents a 
confinement mechanism of the spinors inside the solitons. For a special submodel there ex- 
ists an equivalence between a U{1) vector conserved current, bilinear in the spinors, and a 
topological currents depending only on the first derivative of some scalars. The equivalence is 
such that the charge density associated to that U{1) current can only exist on those regions 
where the space derivative of the scalars are non-vanishing. If one then consider excitations 
of the fields around a given soliton solution, one observes the charge density can exists only 
inside the solitons, since the scalars tend to constants outside them. Therefore outside the 
solitons, the spinors can live in "white" states carrying no U{1) charge. That resembles very 
much what one expexts to happen with the confinement mechanism of QCD, even though 
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our model presents a bag model like confinement instead of the dual Meissner effect. 

The model possesses a zero curvature representation based on the s[3(C) affine kac-Moody 
algebra. It constitutes a particular example of the so-called affine Toda models coupled to 
matter fields which have been introduced in paper The corresponding model associated 

to s[2(C) has been studied in paper where it was shown, using bosonization methods, that 
the equivalence between the currents holds true at the quantum level and so the confinement 
mechanism does take place in the quantum theory. 

The objective of the present paper is to calculate exactly the three species of one-soliton 
solutions as well as the corrensponding six two-soliton solutions describing the scattering 
between any two of them. We also evaluate the time-delays for the interactions among the 
solitons. The character of the interactions, i.e. the relative strength and if it is attractive 
or repulsive, is determined by the scalar product of the roots associated to the solution. 
The solutions are constructed using a mixture of the dressing P] and Hirota methods. 
The dressing method in association to the vertex operator representations of the Kac-Moody 
algebra is quite good in finding the spectrum of solitons, to construct the so-called tau- 
functions 0, and to provide an ansatz for the solution. The Hirota method, which does not 
provide us any hint of the correct tau-functions, is then used to actually evaluate the explicit 
expression of the soliton solutions. We have implemented an algorithm in Mathematica to 
execute the procedures of the Hirota method. In this way we escape from two difficulties of 
the methods. The Hirota method needs the tau-functions and an ansatz for them, but it does 
not tell us how to construct those. That is solved through the dressing method, which in its 
turn needs the evaluation of matrix elements in the vertex operator representation. We avoid 
that calculation through the Hirota expansion method. 

The paper is organized as follows. In Section |^ we describe in detail the model, its 
symmetries and the zero curvature representation. In Section ^ we discuss the dressing and 
Hirota methods, construct the tau-functions and the corresponding Hirota's equations. In 
Section ^ we discuss how to use the infinetely degenerate vacua to introduce non-trivial 
topological charges. Although we work with a slsi^C) Kac-Moody algebra, the topological 

charges take values in the weight lattices of three 5 [2(C) subalgebras. In Section || we show 
how to use the breaking of the conformal symmetry to evaluate the masses of the solitons 
just from their asymptotic behaviour. In Section ^ we show how to evaluate the time delays 
associated to the interactions of the solitons in a two-soliton solution. In Sections |^ and |^ the 
one-soliton and two-soliton solution respectively are calculated. The appendix ^ contains the 
results about Kac-Mooody needed for the calulations performed in the paper, and in appendix 
we present the algorithm for MATHEMATICA used to evaluate the soliton solutions. 

As a matter of conventions, we use for the space-time coordinates in the two-dimensional 
Minkowski space the notations = t and = x. The light front coordiantes are defined 
as = x° ± x^. Hence, we have d± = | (Sq ± (9i), and d'^ = d'^ — df = 49+ (9_. For the 
7-matrices we use the following representation: 

7o = -i(j "J), 7i = -i(i J), 75 = 7o7i=(j _J ) • 
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2 Description of the model 



The affine Toda models coupled to matter fields have been introduced in paper p. In the 

present paper we consider the model associated to the affine Kac-Mody algebra g = slsi^C), 
which equations of motion are[| 

:i+75) .3. (1 -75 )' 
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(1 + 75) 



(1 - 75) 



d^T] = 0, 
iYd^tp' = rUiVi 



(I-75) 



1 ^ 



i=l 



i=l 



(I+75) , 3 JI-75) 

2 ^ 2 

3. (1 + 75) (I-75) 
2 2 



a = 1,2, 



z = l,2,3. 



i = 1,2,3, 



(2.1) 

(2.2) 

(2.3) 
(2.4) 

(2.5) 



where ms = mi +m2, V'* and ip\ i = 1,2, 3, are independent Dirac spinors, and = (■?/'*)^7o. 
Notice that we can not take ip^ to be the complex conjugate of ip^* , since that is incompatible 
with the equations of motion. In addition, we have 



75 (2(^1 -(^2) „75'? 



^i75(2<^2-¥'i) „75''? 



— gi75(¥'l+¥'2) g275??^ 



(2.6) 



Denoting the two component spinors as 



(2.7) 



and similarly for ijj^ and W, we see that 

u}, = ul = ul = 0, 



ul = ul = ul = 0, 



(2.^ 



and 



u 



u 



m2m3 

mi 

mimg 

e ' 

m2 
mim2 



3 Wil p-«(</'i+'^2) 



2 „j(2(^2-¥'i) 



2 A{2^i-^2)] 



m2m3 

e' 

mi 



^2 ^3^.(2^2-^.) +^2 ^3 g 



2 ^-*(</'i+¥'2) 



mims 
m2 

17117712 



^2 ^.(2^2-^1) 



1 ^72 i(2(/3i-v?2) 



(2.9) 
(2.10) 
(2.11) 
(2.12) 
(2.13) 
(2.14) 



^The spinors have been rescaled by a factor iy^mT/i with respect to those in and the scalars cpa by a 
factor i. 
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The model under consideration admits a representation in terms of the Lax-Zakharov-Shabat 
zero curvature condition 



d+A_ - d^A+ + [A+ , A_] = 0, 
where the potentials are given by 

A.=-BF+B-K A_ = -d-B B-^ + P- . 



(2.15) 



(2.16) 



An important ingredient in the definition of the potentials is the use of the principal gradation 
of g (see the definition in the appendix 0). The four scalar fields of the theory live in the 
subgroup B obtained by exponentiating the zero grade subalgebra go given in ( [A.24| ), and we 
parametrize it as 

]^ = (^l=i^faH°+riQpp^i+uC ^ (2.17) 

The model possesses six Dirac spinor fields V'*, and V'*, "^ = 1,2, 3, which live in the subspaces 
with non-zero grades, i.e. we have 



(2.18) 



where i?±3 are fixed elements of the subspaces g±3, and the mappings and F2 take values 
in the subspaces g±i and g±2 respectively. We choose the elements E±2, as 



^±3 = g [(2mi + ma) Hf^ + (2m2 + mi) H^^] 

and use the following parametrisation for and F2'- 
„^ 1 



2 
1 

2 
1 

2 
1 



^3 + ^1 V'lj + m2 El 



02 



mi El + m, El + ms E\^ 
ms i^l E-l - mi V^i E\^ - ma V^i 
mi + m2 V-l - rns V^l ^° , 



(2.19) 

(2.20) 
(2.21) 
(2.22) 
(2.23) 



The model described by equations (|2.1|) - (|2.5| ) is invariant under the conformal transfor- 
mations 

x+^f{x+), x.^g{x.), (2.24) 

with the fields transforming as 



^Pa[,X+ , X 

e 



!fa{X+ , X. 

e~ 



(fia[X+, X-), 



-r){x+ 



if')^^^ {g')^^^ e~^''^+ 



^\x+ , X.) ^ ij\x+ , x_) = e5(i+^5)iog(/')^-^ ^hi^-^^)^osig'r^ 

iIj\x+ ,X^)^ij {X+ , X_) = e^(l+^5)l°s(/')"^ Qh(^-ls)los{g')^-' ^ 



(2.25) 
(2.26) 
(2.27) 

(2.28) 

(2.29) 
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where h (aj) is the height of the root, i.e. h (ai) = h (oj) = 1 and h (a^) = 2, and where 
the conformal weight 6, associated to e"*^, is arbitrary. The two dimensional Lorentz trans- 
formations Xj^ Xx^ and x^ x^/X are contained in ( |2.24| ), as one observes by choosing 
f(x+) = Ax+ and g{x^) = x^/X. 

Since -£^±3 commutes with the Cartan subalgebra generators Hi and if 2? it follows that 
equations (p.l|) - (p.5|) are also invariant under the transformations 

X") hL{x') x') hR{x+), (2.30) 

F+{x+,x-) ^ (2.31) 

F-{x^,x-)^hLix-)F-ix^,x-)hl\x-), (2.32) 

where 

h^ix-) = e'[«i(--)^^+«i(--)^^l, hn{x+) = e'[«H(-^)^^+«i(--')^^]. (2.33) 
That implies that the fields transform as 

^a^^a + ^L + Cn, V^V, (2.34) 

^1 ^ e^i'^+-/r^)i'-^eR+eH) e5(i-75)(2€i-5i) ^1^ (2.35) 

^2 ^ (l+^5)(4_24) g i +2ei) ^2^ (2.36) 

^3 ^ ei(i-^5)«i+5i) ^3 (2.37) 

and ijj^ transforms in the same way as with C,fiL replaced by —C,fiL- 

Notice that, by taking ^% = = with constant e^", we get a GLi(C) x GLi(C) 

global symmetry, with the fields transforming as 

V^a^V^a, V^V, "u^u (2.38) 

and 

^1 ^ e'(2^'-^'V\ V'' ^ e'(2'^'-^'V', V^' ^ e'(^'+^'V' (2.39) 

with i/)* transforming in the same way as ip"^ with 6°" replaced by —6"". 

On the other hand, if we take = ^2 = ~\x°'i with constant x", we get a GLi(C) x 
GLi(C) chiral symmetry, 

</5a V^a - X", z7 (2.40) 

and 

^1 ^ e^75(2x^-x^y ^ ^2 ^ e'^^^^x^-'^'y , ^3 ^ gi75(x^+x^y (2.41) 

with ip^ transforming in the same way as ip"^ with replaced by — x". 

If our theory had a Lagrangian it would follow from the Noether theorem and the above 
GLi(C) X GLi(C) symmetries that it would possess two vector and two chiral conserved 
currents. However, a careful analysis of the equations of motion (|2.1| )- (|2.5| ) reveals that only 
half of such currents exist. Indeed, one can check that, as a consequence of (|2.1|)-(P3|), the 
vector current 

3 

J, = iY,^1i''l,i^' (2.42) 
1=1 
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and the chiral current 

3 

= i ^ ^Sm75^' + 2 df,{mi^i + 1112^2) (2.43) 
1=1 

are conserved 

^^'J^ = 0, 9^JJ = 0. (2.44) 
The existence of these two conserved currents imphes that the quantities 

3 

J = i ^mji/^i/j^ + 2d+{mi<^i + m2<^2) (2.45) 

and 



i=l 



J = \^m\i)\i)\-2d-{m^^^^m2ip2) (2.46) 
1=1 

satisfy 

d_J = 0, d+J = 0. (2.47) 

Indeed, the currents = are obviously conserved. Their hght front components 

have the form 

3 

4+^ = 2 i ^ + 2 5+ (miypi + m2^2) , J-^ = 2d_ {m,^, + 11x2^2) , (2.48) 

i=\ 

3 

= -2 5+(miV9i+m2(^2), /"^ = 2i ^m2V^iV'l-29_(mi^i + m2<^2). (2.49) 



j=i 

and the conservation laws 

+9.4^^ =0 (2.50) 

lead to relations (|2.47| ). 

One can now perform a reduction of the theory by imposing the constraints 

J = 0, ^ = 0, (2.51) 

which are equivalent to0 

3 

i ^ m\ 7/^'7V = -2 e^"a, (mi(/Pi + m2y^2) (2.52) 

Consequently, in the submodel defined by ( p.51|) , the vector current (|2.42|) is proportional to 
a topological current. The relation (|2.52|) has very important consequences for the physical 
properties of the theory. The time component of it implies that the charge density associated 
to the current ( |2.42| ) is proportional to the space derivative of the 93 fields, i.e. 

3 

jO = i ^ ^'-ff = -2 d,{mi^i + 7712^2). (2.53) 



^We use the normalisation e'^^ = 1. 
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Consequently, there can exist charges only on those regions where the space derivative of the 
scalar fields are non-vanishing. As we will see below the soliton solutions of this theory are 
localized in space, in the sense that the scalar fields are not constant only in a region with a size 
determined by the soliton masses. Therefore, in the quantum theory, if we look for excitations 
around a soliton solution, the spinor field can live outside the solitons in configurations with 
charge zero only. We then have a confining mechanism which resembles that of the bag model 
for QCD. 



3 Dressing transformations and the Hirota method 

We construct the soliton solutions using a combination of the so called dressing transfomantion 
method and the Hirota method 0. Let us start with the discussions of the former one. 

The fact that the potentials A± in ( p.l6| ) satisfy the zero curvature ( p. 151) implies that 
they can be represented as 

A^ = -d^TT-\ (3.1) 

where T is a group element obtained by exponentiating the affine Kac-Moody algebra s 13(C) 
(see appendix ^). In addition, the potentials A± satisfy the so called grading condition 
of the Leznov-Saveliev method ||^. This condition states that the potentials A± has some 
limitations on their decomposition over grading subspaces. In the case under consideration 
the potential A_ has only components belonging to the subspaces go, fl-i, 0-2 and 0_3, and 
the potential A^ has only components from the subspaces g+i, 0+2 and 0+3. It can be shown 
@] that potentials satisfying the grading condition can be represented in form ( p. 16 ). Here 
E±3 are abitrary mappings taking values in g±s and satisfying the conditions 

d+E_, = 0, 9_E+3 = 0. (3.2) 



Choosing them in form ( p.l9| ) one fixes the class of equations under consideration. The dress- 



ing transformation method allows to construct gauge transformations which do not violate 
the grading condition and the form of -^±3. Therefore one gets new solutions to the equations 
of motion starting from some initial solution to these equations. 

The concrete procedure to obtain desired gauge transformation is as follows. Consider a 
constant group elements p such that the element 

S = TpT-^ (3.3) 

admits the generalized Gauss decomposition 

E = E_EoS+, (3.4) 

where Uq, and take values in the subgroups corresponding respectively to the subal- 
gebras g<o, 0o and 0>o defined in appendix 0. If decomposition (|3.4| ) exists, we can introduce 
a new group element related to T in two different ways: 

TP = E+T = EZ^T p. (3.5) 

We now define the fiat connection 

AP = -d,TP TP~^ = e^A^O-^^ - 8,0^0^^ (3.6) 
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with ^ 

0+ = r+, = E^^ EZ^ = 0^. (3.7) 

We see that and are related by two gauge tranformations; one with a group element 
involving only non-negative grade generators and the other only non-positive grade generators. 
That fact implies that satisfies the grading condition. Moreover, it can be shown that 
the gauge transformations under consideration preserve the form of E±^. Consequently, the 
dressing method defines transformations on the space of solutions of the model defined by 
the zero curvature condition ( |2.15| ). Such transformations are parametrised by the constant 



group element p. 

A quite general procedure to construct soliton solutions in integrable theories, using the 
dressing transformation method, is described in paper It constitutes a generalization to 
practically any two-dimensional integrable hierarchy, of the so-called "solitonic specialization" 
in the context of the Leznov-Saveliev solution for Toda type models |P, 0, . The idea is to 
start from a "vacuum" solution such that the connections A±, when evaluated on it, belong to 
an algebra of oscillators, i. e. an abelian (up to central terms) subalgebra of the Kac-Moody 
algebra. One then looks for the eigenvectors V^, in g, of such oscillators. The solitons belong 
to the orbits of solutions obtained by the dressing transformation performed by elements of 
the form p = e^'i e^^a . . . e^*" . 

Notice that 

1 ^ 

-— ^m^x+x", (^„ = 7^ = ^* = V^^ = 0, a = 1,2, i = 1,2,3, (3.8) 

i=l 



is a "vacuum" solution of the equations of motion (|2.1|)-(|275|). The potentials ( |2.16| ) evaluated 
on it become 

1 ^ 

/l+ = -^+3, A_ = E_, + — J2'^'^^^ (3.9) 

1=1 

and so we can take the initial group element T in the forniQ 

T = e^^^+3 e-^"^-^ (3.10) 



Here we have used the equality 



[^+3,i^^-3] = ^X^m^C (3.11) 



12 

i=l 



which follows from (|2.19|) and ( |A.12| ). 



Notice that E±3 belong to the "algebra of oscillators" generated by the elements (see 
relation ( |A.12| )). That is the so-called homogeneous Heisenberg subalgebra of the affine 

Kac-Moody algebra sl^i^C). 
The eigenvectors of £^±3 are 

00 

V±^Az)= J2 ^ = 1,2,3, (3.12) 



^Recall that the group element T entering (|3.l| ) is defined up to the right multipHcation by a constant 
element of the group. 
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with z being an arbitrary parameter. Indeed, one has 

y±.Xz)\ = ±\my^^X^), [i?_3, y±.Xz)\ = ±^m,V^±,,(2;). (3.13) 

Notice that VaX^) cind y_Q,.(— 2;) have the same eigenvalues. It turns out that the soliton 
solutions are obtained by taking the constant group element p to be products of exponentials 
of the operators 

V^X<. ^) = <V^A^) + a^y_^X-z). (3.14) 



We now perform the dressing transformation (|3.6| ) starting from the vacuum potential 



]9|). The transformed potential has the structure given by relation ( 2.16 ) with B and 
given by ( p.l7| ) and (|2.18|) , therefore we have the equality 

b{E+s + + F^)b-' = 0±E+sOg' + 9+6»± Og\ (3.15) 

where 

5 = gi(VlHl+^2H2)^ (3.16) 

and the equality 

- i 5_(^i Hi + ^2 H2) -d_uC + E_3 + + 

1 ^ 

= O^E^^Og' - 0^' + -J2 C. (3.17) 



i=l 



These equalities relate the new fields and the parameters of the dressing transformation. 
Notice that the rj field appears in ( p.l6|) multiplying the grading operator Qppai, and that 
involves the operator D (see ( A. 201) ). Since D is not the result of any commutator it follows 



that the dressing transformation method does not excite the field r] if one starts from a 
solution where it vanishes. 

In order to construct the solution solutions we have to split ( p. 15] ) and ( p.l7| ) into the 



grading subspaces. For instance, taking grade 3 component of (|3.15| ) for and the grade 
one of ( |3.17| ) for 0_ too, gives 

Tq-^ = 66(^+^^-1"^'"^"")^. (3.18) 
Then from ( p.4|) and ( p.lOj ) one gets for 



The idea now is to consider matrix elements of both sides of ( p.l9| ) in states of the funda- 



mental representations of q to obtain the explict space-time dependence of the fields evaluated 
on the solutions determined by the constant group element p. For instance, one obtains from 
(|A|) that 

To To 

with 

r, = (A,|r|A,), J = 0,1, 2, (3.21) 
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where | Xj ) are the highest weight states of the fundamental representations of g, satisfying 
( |A.28| ) and (|A.29|) , and where we have used that, as a consequence of ( |A.29| ) and (|A.33|) , 



0+|A,) = |A,), {X,\0 



-1 



(A, 



(3.22) 



In order to obtain the expressions for the spinor fields we represent the mappings 0+ and 

0~ as 

= eS™>o , o_ = eE™>o , (3.23) 
where the mappings and t^^™) take values in the subsapces 0+^ and Q-m, respectively. 



According to ( |A.29| ) and ( |A.33| ) for the highest weight vectors of the fundamental reparesen- 
tations we have 



I A^. ) = 0, ( A, 1 = 0, 



m > 0, 



J 



0,1,2. 



(3.24) 



Considering the components of grades 1 and 2 of (|3.15D for implies that 



^2 



1 



Analogously, considering the components of grades 

F^ = [t^'\ E^,] ; Ff = [t^'\ E_3] + ^[t^'\ [t^'\ E 



(3.25) 



1 and —2 of (|3.17|) for 6*+, one gets 
1 



(3.26) 



Now, multiplying both sides of (|3.19|) from the left by generators of grades 1 and 2 and 
taking the expectation value on the states \ Xj), one can relate the components of t^~^^ and 
t*-"^-* to the matrix elements of E. The peculiarity here is that to extract t^~'^^ one needs to 
consider two matrix elements for each one of its components. By adding and subtracting the 
corresponding relations coming from (|3.19|) , one obtains the components of in terms of 
matrix elements of E, and a quadratic relation involving the components of t^~^\ A similar 
thing happens when we multiply both sides of (|3.19|) from the right by generators of grades —1 
and —2 to relate t^~^^^ and t*-'''^^ to matrix elements of S. Once we have obtained the relations 
between and t^^^) matrix elements of E, we can use (|]2|), ( ^lEj ) and (F^-(^) 
to obtain the explict space-time dependence of spinor fields on the solutions determined by 
p. The results are 



i'R 



4 



^R,(2) 


m2 ^R,{0) 










(3.27) 


mi T2 


mi To 


Tl 










^1 = 








(3.28) 


1712 n 


m2 To 


T2 






^R 




^1 = 


m2 ^L,(l) 
mg Tl 


mi 
+ — 

mg 


T-"3 
^L,i2) 


(3.29) 


^R 




= 


m2 ^L,(0) 


mg 




(3.30) 






mi To 


nil 




^R 




^1 = 


mi ^L,(o) 


ms 




(3.31) 


r2 




m2 To 


m2 


Tl 


as 

m2 ^R,(l) 


mi T^R^(2) 


^1 = 


~Q!3 






(3.32) 




ms T2 


^0 ' 
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where tj, j = 0,1,2, are defined in ( p.21|) , and where we have denoted 





= (A2I 


^-oi^l A2), 




= (Aol 


^ioi^l Ao), 


(3.33) 




= (All 


^^°aj Al), 


R,(") 


= (Aol 


E^^^S Ao ), 


(3.34) 




= (All 


^-02^1^1)' 


T-0!2 

' T 

l_i 


= (A2 


^^-02 1^2), 


(3.35) 


It 


= (Aol 


E Ao ), 


' T 10\ 


= -(A 


2 UE^_ A2 ), 


(3.36) 




= -(A 


i\EEWX,), 


~ai 

iX 


= (Ai 


^ai^lAl), 


(3.37) 




= (A2I 






= (Ao 


^^ai'l Ao), 


(3.38) 




= (A2I 


Ea,^\ \ A2), 




= (All 


^^^,MAi), 


(3.39) 


^L,(0) 


= (Ao| 


^^."2'IAO), 


^"3 

^H,(2) 


= (A2 


<^|A2), 


(3.40) 




= (Ai| 


<^|Ai), 


0,3 


= (Ao 




(3.41) 



As we mentioned above in the process of obtaining the components of t*^^^^ in terms of the 
matrix elements of E one gets quadratic relations involving the components of t'^^^\ These 
relations are given by 



~Q,2 


^R 


+ ^2r-RA.) 


-^oT^;(2) 


= 0, 


^R 


^R 


'1 ' RAO) 




= 0, 


(3.42) 


Tr' 


^"2 
^R 


-^^^RAD 




= 0, 






~ ^0 ^L,(2) 


+ ^2 ^lIo) 


= 0, 


(3.43) 




rT 


+ ^o^r(i) 


- ^1 ^LAO) 


= 0, 




a2 




- ^ ^1(2) 


= 0. 


(3.44) 



According to the "solitonic specialization" 0, [ly, ^ , the soliton solutions are obtained 
by taking p to be the product of exponentials of the operators (|3.14| ). The one-soliton solutions 
are obtained by taking 



p = p. = exp[\4,(a„^,2;)], z = 1,2,3. 



(3.45) 



Therefore, there are three species of solitons. Similarly, the two-soliton solutions are obtained 
by taking 

P = pi Pj = exp [Va, {a^^ , zi)] exp [14^. {b^. ,Z2)], (3.46) 

and so there are six of such type of solutions. Using ( |3.13|) one then gets from ( p.l9|) that for 
the one-soliton solutions 



E = Ei = e^''^+3 e"^ ^- 
and for the two-soliton solutions 



exp [e^'(^)\4,(a^ 



(3.47) 



Si Si 



where 



with 



7 



Ei iz] 



exp [e^'(^i)v;,(aj,^0] exp [e^^^^-^V^^{h%, z^) 

= 7 (x — vt) 



■ rrii 



rrii 



X 



zx — 



(sign z) 



errii cosh 6' ; 



^2 + 1 



tanh^, 



(3.48) 
(3.49) 

(3.50) 
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where we have introduced the rapidity 6' as z = ee^, with e = ibl. Therefore, if z is real we 
have |f I < 1, where the speed of hght has been normahzed to unity. A generahzation to the 
n-sohton case is evident. 

In order to obtain the final expressions for one-soliton and two-soliton solutions one has 
to evaluate the matrix elements in ( |3.21| ) and ( |3.33| )-( |3.41| ) for the group elements ( |3.47| ) and 
( p.48| ). That calculation can be performed using vertex operator representations for g [|^, ^ 



We have to use in fact the homogeneous vertex operator realization of the three fundamental 
representations of g. It turns out that such representations are integrable in the sense that 
the step operators of g are nilpotent [rT|. Indeed, one can verify that the operators ( 3.12 ) 
satisfy [|12], ^ 

(^i) Kv. (^2) -^0 as 2:1 ^ Z2. (3.51) 

That implies that the matrix elements in ( p.21| ) and (|3.33|) -( ^.41|) truncate when one expands 
the exponentials in (p.47|) and (|3.48|) , and consequently they become polinomials in the c^^'s. 
In this sense, those matrix elements are the Hirota tau-functions for the model described by 
equations ( |2.1| )-( |2l5| ). Note that due to ( p.51| ) the maximal power of each e^' in the expansion 



of a tau-function over them is equal to 2. Therefore, a general tau-function for an n-soliton 
solution has the form 



r 



E 



Plv,Pn = 



-Pl.-.p, 



exp [piri,(zi) H Wn^iS^^X 



(3.52) 



The connection between dressing transformations, solitonic specialization and the Hirota 
method has been discussed in for any hierarchy of integrable models possessing a zero 
curvature representation in terms of an afiine Kac-Moody algebra. It is a quite general result 
that explains a lot of the structures known in soliton theory. In addition, on the practical 
side it allows one to use the good features of each method. 

The Hirota method is quite powerful to actually evaluate the explicit expression for the 
soliton solutions. Besides, being a recursive method it allows a simple implementation on 
a computer program for algebraic manipulation like Mathematica. However, one of its 
difficulty is to find the relation between the tau-functions and the fields that lead to the 
truncation of the Hirota expansion. But the dressing transformation method together with 
the solitonic specialization does exactly that. 

The dressing method on the other hand, requires the evaluation of matrix elements in 
vertex operator representations of the Kac-Moody algebra, which in the case of higher soliton 
solutions become extremely tedious. One can take advantage of those features to speed up 
calculations. 

To use the Hirota method we have to find the equations satisfied by tau-functions. They 
are determined by substituting the relations between the fields and tau-functions ( p.20|) and 
( p.27] )- (|3.32| ) into the equations of motion (|2.1| )- (|2.5| ). Notice that, except for 77, we have 15 
fields in the model ( p.l| )-( pl5[ ) (three scalars and six two-component spinors). However, in 
( p.20| ) and (|3.27|) - (|3.32|) we have introduced 21 tau-functions. The six equations missing are 
the algebraic relations ( p.42|) -( ^.44| ) provided by the dressing method. 

Substituting ( p.20|) and (|3.27| )- (|3.32|) into ( |2.1|) -( ^75D one gets that the equations for 
and V lead to the following three Hirota equations 



4 Ti T2 (ro 9+9_ro - (9_ro S+tq) 



12 



- mi m2 Ti Ta r°^(,) f^^(2) - mi^ r^^^) ^^2) ' (3.53) 
4 To r2(ri 9+9_ri - 9_ri S+n) = ma' ri f^^ + mg' tq Ti r^^ 

- m2 ma Tq t'^^^^^ + m^" r^^ r^;^^) 

+ ma' To' f^;(2) ^^,'(2) - "^2 To r2 f^;(o) r°^(2) , (3.54) 
4 To ri (r2 9+9_r2 - 9_r2 9+r2) = mi' n r2 f^^ + mg' tq r2 r^^ 

- mi ma Tq n f^^^^^ r^;^^^) + mi' n' f^^^^^ t^^^^^ 

+ ma' To' r^_^(i) - mi ma tq Ti f^^g) r^^^^.^ . (3.55) 

The components of the equations of motion (|2.1| )- (|2.5| ) for -i/^'s and t/^'s that do not involve 



the quadratic terms f/* and V (see ( p.8| )) lead to the Hirota equations 

2 (n a+rf - 9+ri) = ma tq r£(2) - m2 r^^^y (3.56) 

2 (r2 S+rf - rf a+r2) = ma To rl\^^ - n t^^^^,^ (3.57) 

2 (ro a_r;^^ - r]^^ a_ro) = m2 r^^,) + mi n r^^^), (3-58) 

2 (n - a_ri) = ma To r^]^,^ - m2 fl]^,^ (3.59) 

2 (r2 - f- a_r2) = ma To f^l,^ - n f-„), (3.60) 

2 (ro S+ff - 9+ro) = -m2 f^;^,^ - m, n f"^^^)' (3-61) 



and the other components of the equations ( p.lD -( ^^ for ip^s and ip^s involving the quadratic 
terms and lead to 

2 mi tI (r2 9+r^|2) - r'^^h + 2 m2 r^' (n S+r^^^^^ - r^^^^^ S+n) 

= "^2 ^2 ^r\q) - ^0 ri T2 - m2 ma tq r2 r^^ r;^i(2) 

- < Tll^^ + mi ma To n r^^^ , (3.62) 

2 m2 r| (ro S^r^;^,) - r^^^^ d^r,) + 2 ma r^' (r^^^^^) ^-^2 - r2 d^r^^^,^) 
= ml To Ti T2 t2' - ml rl f^^^^ + mi ma tq Ti r'^ f "^(g^ 

- ^2 rl ft ^£(1) - "^1 ^2 Ti r2 f^^ Tll^^ , (3.63) 

2 mi rl (ro a_r-,) - r-,) a_ro) - 2 ma To' (n 9_r- - r- S^n) 
= m\ To Ti r2 + m^ To' r^^ - ma To t'^ f^^^^ 
+ mi m2 ri r2 f^^ r^"^(,) + m? Ti' r^^^) , (3-64) 

2 m2 r2' (ro d+f^^^^ - f^]^^^ S+tq) + ma {f^^^^^^ d+T2 - T2 d+f^h)) 
= ml To Ti T2 + mi ma Tq n t^^^^ - mj r^'^^^^ 
+ "^2 ^2 ^Si) + ^1 ri T2 rl- f2{2) ' (3-65) 

2 mi Ti' (ro - f^lo) ^+^0) + ma (f^^^^.^ S+n - n 9+f£|^) 

= m^ To Ti T2 - m2 ma Tq r2 f °« r^_\o) + mj t'^' r^;(2) 

- mi m2 ri r2 rf f^iX) ' ^^2) ' (3-66) 
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= ml To Ti Ta f^^ + ma mg Tq f £^(2) " "^2 ^1 ^^'(0) 

- mi ms ro n + m? r^^ . (3.67) 

The equivalence between the vector and topological currents (|2.52| ) in terms of the tau- 
functions become 

2 mi T2 (n 9+ To - To d+Ti) + 2 ma Ti (rg O+Tq - Tq S+ra) 

= mi m2 T2 r^;(o) - mi Tq r^i(2) ^fl' + ^2 n r^_'(o) 

- ma Tq r^J'^^) r]^' - ma mg ra r]^' f]^_'(^ - mi mg n r]^'' f^|'(2) ' (3-68) 
2 mi ra (n 9_ro - tq d^n) + 2 ma n (ra 9_ro - ro 9_ra) 

= ma ma ra r^^ f^^ + mi mg ri r£=|2) ^2' + ^1 ^2 ^2 ^£'(0) 

- mi mg ro r^^ f^^ 2) + "^1 "^2 n ^£'0) ~ "^2 "^3 ^0 1"^' ^L,a) ' ^^'^^^ 

Resuming we can say that the soliton solution can then be obtained by evaluating the 
matrix elements ( p.21| ) and (|3.33| )- (|3.41|) in the homogeneous vertex operator realization of 

the three fundamental representations of the affine Kac-Moody algebra s[3(C). Alternatively 
one can obtain the same solutions by using the Hirota method. In such case one uses the 
ansatz ( p.52| ) for the tau-function provided by the dressing method. The coefficients Cp^...p„ 
can be determined recursively using a program for algebraic manipulation. In appendix ^ we 
give the code of a Mathematica program implementing the Hirota method. The explicit 
soliton solutions are given in Sections |^ and ||. 



4 Topological charge 

In (1 + 1) dimensions the topological current is defined as J*°p ~ e^^^d'^^, where $ is some 
Lorentz scalar quantity. The reason is that if ^ does not present singularities, the space- 
time derivatives acting on it commute and therefore J*°p is conserved independently of the 
equations of motion, i.e. 9^ J*°p = 0. The topological charge is then = J dx Jq°'^ ~ 
${+oo) — <P{—oo). For finite energy solutions the fields have to approach a vacuum con- 
figuration for X ±00, and therefore the introduction of a topological current only makes 
sense in theories with degenerate vacua. Our theory ( |2.1| )-( pl5| ) has four scalars fields from 
which we can build topological charges. However, rj is a free field and its vacuum configura- 
tion do not lead to interesting charges. The field u is not suitable either because one needs 
u — > —j^ Ylt=i as X — s> ±00. Therefore, we are left with the two scalars (pa, a = 1,2. 

Then, following the approach adopted in the abelian affine Toda models we introduce the 
quantity 

f-t'-^f-. (4.1) 

a=l « 

where a^, a = 1,2, are the simple roots of sls{C) (see appendix We then have that 



(Pa = {(p\Xa), where are the fundamental weights of s[3(C) defined by the relation |jT3 



2^ = 5... (4.2) 
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One observes from (|2.9| )- (|2.6| ) that the fields enter into the equations of motion through 
the combinations 6^*^*^'"^^ j = 1, 2, 3, which are invariant under the transformations 



V9 — i> (y9 + 27r/i 



(4.3) 



where [i is any vector on the root space satisfying the condition (ai|/i) G Z. Such vectors 
are the so-called weights of s[3(C) and they constitute an infinite discrete lattice called the 
weight lattice |]13i. In the case of the abelian affine Toda models the vacua of the theory is 



determined by such weight lattice. However, in the model described by equations (|2.1| )- (|2.5| ) 
one notices that for any constants yji"^ and Vj^^^ 



0, 



'-Pa 



1 ^ 

-y 

12 ^ 



m, X X 



(4.4) 



is a vacuum configuration. Consequently, the vacuum configurations of the (pa fields are not 
determined by the weight lattice of s[3(C). What we will see below is that the topological 
charges of the one-soliton solutions of (|2.1| )- (|2.5| ) lie in fact on the weight lattice of the three 
s[2(C) subalgebras associated to the three positive roots of 513(C). 
We shall define the topological current and charge as 



27r '^'^ 



d'ip, Q 



top 



dx Jo°^ 



27T 



(4.5) 



In terms of the tau-functions introduced in ( |3.20|) one can express the topological charge as 



Q 



top 



1 

2^ 



+00 



(4.6) 



The relation ( p.20|) implies that —if a = ln|Ta/ro| + iarg(ra/ro). Therefore, in order to have 
real solutions for the fields ipa, one needs \Ta\ = \tq\. Given the solutions for the tau-functions, 
the relation (|3.2CI|) determines ipa only modulo 2tt. Therefore, the topological charges are 
defined up to an integer linear combination of the simple roots aa (recal that we use the 
normalization = 2). When evaluating the charges we should then use the prescription 
that the arguments of the tau-functions lie between vr and — vr. 

Notice that the topological current in (|2.52|), which is equivalent to the spinor vector 



current is the projection of ([4.5|) onto the vector —4tt {rriiXi + 7712X2)- 



5 Masses of fundamental particles and solitons 



The model described by equations (|2.1| )- (|2.5| ) is conformally invariant. It means, in particular. 



that if one has a solution of the equations with some finite mass, then by a scale transformation 
it can be continuosly transformed to a solution with any other mass. Hence to have stable 
solutions we have to break the conformal invariance. The simplest way to do this is to freeze 
the field rj. This fixes also the masses of the fundamental fields. Indeed, taking, for example, 
77 = 0, and considering the linear part of the equations ( p.l|) -( ^3D one observes that the 
masses of the fundamental particles are (see section 7 of 0) 

m^i=m^i=mi, m^^=m^ = Q. (5.1) 
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The mass of a soliton is defined to be its rest frame energy. The energy is measured 
by the Hamiltonian or the energy-momentum tensor. The model under consideration does 
not possess a Lagrangian and so we cannot calculate the canonical energy-momentum tensor. 
However, in reference it was shown that our model belongs to a class of theories obtained by 
the Hamiltonian reduction from the two-loop WZNW model ||T^ |T5[. The energy-momentum 
tensor L^,^, of the two-loop WZNW model is of the Sugawara form, i.e. quadratic in the 
conserved currents. We take the energy momentum tensor of our model to be the reduced 
form of Lfj^i, under the Hamiltonian reduction mentioned above. Denote that by Uuy- Such 



tensor is conserved and traceless, i.e. 



red 



0, 



' red A* 



0. 



(5.2) 



However, according to the arguments given in |T^, § the masses measured by such tensor 
must vanish, since it corresponds to a conformally invariant field theory. On the other hand, 
the conformal symmetry can be broken by freezing the rj field to a constant value. We define 
the energy momentum tensor of such theory with a broken conformal symmetry by 







broken 



= o 



where 
with 



e 



' red 



+ 5, 



(5.3) 



(5.4) 



(5.5) 



and K being the coupling constant of the two- loop WZNW model g |T4|, . The tensor 

is conserved but not traceless. We define the soliton masses to be the space integral of 6>qq°'^'^'^ 

on the soliton rest frame. 

Therefore, using (|5.4|) and (|5.5|) , we obtain the following expression for the mass M of the 
soliton moving with the velocity v. 



M 



dx (6», 



broken 
00 



_ 1 , , X^X 

+ 3 (V^i + V2) + 



E 



j=0 



+00 



(5.6) 



Here we have used relation ( p.20| ), and the fact that the integral of U^^^- vanishes, as dis- 
cussed above, due to the conformal symmetry. The term is due to the nonzero vac- 
uum energy density, associated to the fact that the vacuum configuration of the field z7 is 
— (1/12) Y^i=i . Thus, the soliton mass is determined by the asymptotic behaviour of 

the solution. 



6 Time delays 

A soliton is a classical solution that travels with a constant speed without dispersion and 
keeps its form under scattering on other solitons. The effect of scattering is only a phase 
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shift or a displacement in the sohton position. Now we will show that the solitons we have 
been working with are true solitons in this sense and calculate the so-called time delays of 



the scattering of two solitons, using the procedure elaborated in papers |T^, |§1 • 

Consider two solitons that in the distant past are well apart, then colide near t = 
and then separate again in the distant future. Therefore, except for the region where the 
scattering occurs, the solitons are free and so travel with constant velocities. Let us write the 
trajectories of one of the solitons before and after the collision as 

X = vt + x{I) and x = vt + x{F), (6.1) 

since the velocity is the same. The lateral displacement at fixed time is measured by 

A{x) = x{F) - x{I), (6.2) 

and the time delay is defined as 

A{t)=t{F)-t{I) = -^ (6.3) 

with the intercepts of the trajectories with the time axis being given by t{F) = —x{F)/v and 
t{I) = —x{I)/v. The lateral displacement and the time delay are not Lorentz invariant, and 
so we consider the invariant 

EA{x) = -pA{t) (6.4) 

with E and p = vE being the energy and momentum of the soliton respectively. Since E is 
positive it follows that A{x) has the same sign in any reference frame, with only its value 
being frame dependent. The time delay on the other hand may change the sign under Lorentz 
transformations. One can show that the lateral displacements and time delays for the two 
solitons participating in the collision have to satisfy |]l^ H 

EiAi{x) + E2A2{x) = 0, PiZ\i(t) + P2A2{t) = (6.5) 

where the index i labels the quantities associated to the z-th particle. Notice therefore that, 
since the energies are positive, the lateral displacements have opposite signs. Clearly, in the 
center of momentum frame where pi+P2 = 0, the time delays are equal, i.e. Z\^™(t) = Z^2™(^)- 
Therefore, from (U) we have that E'^'^ Af^ (x) / pi = -E^'^Al'^{x)/pi = -Af^{t) = -Z\r(4- 
Consequently, since A{x) has the same sign in any frame, it follows that, if the first particle 
moves to the right faster then the second particle (so that pi is positive in the center of 
momentum frame), then {—Ai{x)), A2{x), A™{t) and Z\2™(t) all have the same sign. The 
physical interpretation of that sign is related to the character (attractive or repulsive) of the 
interaction forces. Indeed, if the force is attractive then the first particle will accelerate as 
it approaches the second particle and then decelerate. That means that Ai{x) is positive 
and so the common sign negative. Therefore, attractive forces lead to a negative time delay 
in the center of momentum frame, and clearly repulsive forces lead to a positive time delay. 
These considerations assume that the two particles pass through each other, and there is no 
refiection. However, when the masses of the two particles are equal there is the possibility of 
occurring refiection. 

In order to evaluate the time delay we choose the first particle to move faster to the right 
than the second particle, i.e. Vi > V2, with vi > 0, and therefore from ( p.50|) \zi\ > \z2\ or 
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9i > 02, where we have denoted Za = eaC^" with ea = ±1. Let us track the first sohton in time 



17|| , i.e. hold x — vit fixed as time varies. One gets 



X — V2t = X —JJit +{Vi — V2)t. 
const 



(6.6) 



We then get from (|3.49|) that, if 62 = 1, e-'"*^^^) as t — > —00, and e^'^^'^'^ — > 00 as t ^ 00 



{r stands for one of the three Fi of ( |3.49|) , corresponding to the second sohton). For the case 
€2 = —1 the hmits get interchanged. 

As we will see in section taking 62 = 1, all the two-soliton solutions become, in the limit 
t —00, a one-soliton solution. Now, in the limit t 00, the two-solitons also become a 
one-soliton but with the replacement 



Zi - iJZ2 
Zi + UJZ2 



(6.7) 



where n is either 1 or 2, and uj = ±1, depending upon the type of two-soliton solution under 
consideration. Therefore, the relevant effect of the scattering on the solutions is a lateral 
displacement of the first soliton given by 



7i (x - fit) 7i 



X 



fit H In 

7i 



Zi - UJZ2 
Zi + UJZ2 



(6.^ 



If one takes €2 = — 1 instead, one observes that the direction of the arrow in (16. 81) reverses. 



Therefore, using (|6.2| ) one sees that the lateral displacement for the first soliton is given by 



AAx) 



£2 

71 



In 



Zi - IJJZ2 
Zi + LUZ2 



eie2 



m cosh 9i 



In 



.i9i-e2)/2 



eie2Uie 



-{ei-e2)/2 



6162000' 



-{ei-e2)/2 



(6.9) 



where m stands for one of the masses m^, i = 1,2, 3, corresponding to the first soliton. Since 
6162UJ = ±1, one observes that Ai(x) can in fact be written as 



u n 



m cosh 9i 



In 



tanh 



9i — 02 



(6.10) 



Using ( |6.3|) and (|6.10|) , one gets that the time delay is given by (assuming fi > f2, with 
vi>0) 



uj n 



In 



tanh 



6i — 6-2 



(6.11) 



m sinh 9i 

Notice that we have taken 61 > 62 and so the hyperbolic tangent can vary from to 1 and 
therefore its logarithm is always negative, and so from ( |6.10| ) one sees that Ai{x) for fi > f2, 
with fi > 0, has the same sign as u. Therefore, from the above considerations we conclude 
that the forces between the solitons are attractive for uj = 1, and repulsive for uj = —1. In 
section R we discuss the details of each type of two-soliton solution. 



7 One-soliton solutions 

As we have discussed in section H the one-soliton solutions are obtained by taking the constant 



group element p, parametrising the dressing transformations, as in ( 3.45 ). Therefore, we have 
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a species of solitons for each positive root of sl^lC). We give below the results obtained by 
either evaluating the matrix elements ( p.21 ) and ( |3.33D -( ^.41 ), or by using the Hirota method 



(see program in appendix 0). As we show below the one-soliton solution associated to a given 
root Oi excite only the tau-functions (or fields) associated to the same root. Therefore, each 
species of solitons belongs to one of the three s[2(C) subalgebras associated to the positive 
roots of s[3(C). 

We have checked that the three one-soliton and six two-soliton solutions constructed in 
this paper satisfy the equivalence between the vector and topological currents given by (|2.52|) , 
or equivalently ( |3.68| ) and ( |3.69| ). Therefore, they are solutions of the submodel defined by the 



constraints ( p.51|) , and that presents the confinement of the spinor charge inside the solitons. 
7.1 One soliton of species oli 

In this case we have E = Si = exp [e''"^*^^-*V^j (a^^, z)] . The only non- vanishing tau-functions 
are given by 

^„ = l_^^e2^^W, ri = l + ^^e2^^(^), r2 = 1 - e^^^W, (7.1) 

^RM = < ^Rf2) = = < e^^^^^ (7.2) 



.... „. a; 

Z 



As discussed below ( [4.6|) , in order to have real one needs \tc\ = \to\, and consequently 
has to be pure imaginary. The topological current ( |4.(j| ) is then 

Q*°P = -^ai(sign^) (7.4) 

and so it is (up to a sign) the fundamental weight of the s[2(C) subalgebra generated by Hi 
and E±ai. Using ( pl6| ) and ( p.50| ), one obtains that the mass of the soliton is M^"^^^ = 2Kmi. 



7.2 One soliton of species 0:2 

In this case we have E = E2 = exp [e^^^^^Va2iO't2^ ^)]- '^^^ non-vanishing tau-functions are 
given by 

^Q = l_^^e2A(.)^ ^^ = l_^^e2A(.)^ ^2 = l + ^^e2^^(^), (7.5) 

' RAO) ~ OL2 ^ ■> 'Rd)~"a2^^ ' 'L ~"a2^ ' V'"/ 



S2.pr2{z) ~a2 _ _^ 



^:Q2 _ - r2(2) ~a2 __!^„r2(2) ~a2 _ pr2(z) 



Similarly as above, in order to have real we need a^^ a^^ be pure imaginary. The topological 
current is then 

Q*°P = -^«2(sign^) (7.8) 
and so it is (up to a sign) the fundamental weight of the s[2(C) subalgebra generated by H2 

f(a2 
'sol 



and E±a2- The mass of the soliton is m!:^-,^^ = 2Km2. 
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7.3 One soliton of species 0:3 

In this case we have E = S3 = exp ^e'"^^^^Va^{a'^^, z)'j. The non-vanishing tau-functions are 
given by 

z 



_ ^- 0^3(2) -as _ „- „r3(2) ^:"3 03 r3(2) l'7 11'l 

^i?,(l)~"a3^ ' ^i?,{2) ~ "03 ^ ' T^ — ^ ^ • 



In order to have (fa real we need a^^ a^^ be pure imaginary. The topological current (|4.6| ) is 
then 

Q*°P = -^«3(sign^) (7.12) 

and so it is (up to a sign) the fundamental weight of the s[2(C) subalgebra generated by 
Hi + H2 and E±q.^. The mass of the soliton is M^^^^^ = 2Kms. 



8 Two-soliton solutions 

We now present the six two-soliton solutions corresponding to the pairs we can form with the 
three species of solitons constructed in section |^. The solutions are constructed by performing 
the dressing transformations with the constant group element p leading to the elements Eij 
given in ( p.48|) . Instead of evaluating the matrix elements in (|3.21|) and ( p.33|) -( p.41| ) using 



the homogeneous vertex operator realization of the three fundamental representations of the 
0, we calculated the solutions using the Hirota method with ansatz (|3.52| ). The evaluation 



of the corresponding coefficients was performed using the Mathematica program described in 
appendix ^. 

8.1 Two solitons of species ai/ai 

Here we give the two-soliton solution obtained by taking the group element E = Eu = 
exp [e''"^^^^^V^^(a^^, 2:1)] exp [e''"^''^^''V^^(6^^, 2:2)] • The non-vanishing tau-functions are 

^ ^_ ^ai^ti 2 A (21) _ KiKi ^2 A (22) _ i'^ai^ai + ^aiKi) ^1^2 A(2i)+A(22) 

° 4 4 {zi + Z2y 

QaiQgi^ai Ki (^1 ~ ^2)^ 2A(2i)+2A(22) (o i) 

16 (21 + Z2y ' ^ ■ ^ 

^ _^ Qgi Qgi ^2 A (21) _^ ^Qi ^2 A (22) _^ ('^ai Ki ^1 + "ai Ki 4) ^Af^Q+Af^a) 
4 4 (Zi + ^2)^ 

^ Qgi Qgi ^Qi ^ai (^1 ~ ^2)"^ ^2 A (2i)+2 A (22) ^ (^g^2) 



^ ^ _ 2A{2i) _ ^ai 2 A (22) _ (^Qi ^ai + ^ai ^Qi ) ^1^2 A(2i)+A(22) 

' 4 4 (^1 + ^2)' 

^ Qgi ^Qi ~ ^2)^ ^2A(2i)+2A(22) ("g 3\ 

16 (Zi + Z2Y 
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Q«i < K (^1 ~ ^2)^ 

\2 



I '^ai ^ai ^ai 



_ 0^(21) , ,+ A(^2) I ^2)^ 2n(z,)+n(z2) 



I '^ai ^01 ^01 

(^1 - ^2)^ ri(2i)+2ri(^2) 10 p:\ 



~ai ^ gAC^i) _ ^ gri(;.2) _^ ^ai (^1 ^2)^ gA(^i)+2 A(^2) 



i,(2) 



42:1 (2:1 + 2:2)2 



_^ Qgi Qgi ^ai (^1 ^2 A {^i)+A (^2) ("3 g^ 

4^2(^1 + ^2)^ 



~ai ^ _^gA(^i) _ ^gA(^2) _^ '^ai ^ai ^ai (^1 ^A(^i)+2A(^2) 

zi 2:2 42:1 (2:1 + 2:2)2 



^L,(0) 



_^ Qgi Qgi ^ai (^1 ^2)^ ^2 A{^i)+A (^2) /g 
42:2(2:1 + 2:2)2 

^"1 _ pA(^i) , /,+ . „A(22) _ Q«i ^1 (^1 ~ ^2)^ A(^i)+2 A(22) 

\R,(0) ~ Qi 1 ^2 « w , \o 



4(^1 + ^2)2 

Qgi ^ai (^1 - ^2)^ ^2 2 A(^i)+A(^2) 

4 (^1 + ^2)2 



rSVo^ = a+ zie^^(^^) + 6+ z^^^^^^^^ _ <.^^L^WfL^^ pA(..)+2 a(.2) 



iJ,(2) --ai-l- --.1-2- 4(;2l + ^2)^ 

^'^^ " '^'^^''^%2A(.i)+A(.2)^ (8.9) 



Q«i (^1 - ^2)^ ^2 2 A(^l)+A(^2) 



4(^i + ;22)2 

One can check that by taking the hmit e^^^^^^^ 0, such solution becomes the one-sohton 
solution (|7.1|) - (|7.3|) . Now, considering the ratios of tau-functions appearing in relations (|3.20| ) 
and ( |3.27| )-( |3.32| ) one observes that in the limit e^'"^*^^^^ — >■ 00, such solution becomes the one- 
soliton solution ( [7.1|) -( [f73D again but with the replacement 



^ l j e'^^"^^ (8.10) 

and with the ratio a/to changing sign. Therefore, comparing with ( |6.7| ), we observe that the 
lateral displacement and time delay for this two-soliton solution is given by ( |6.10|) and (|6.11|) , 
respectively, with n = 2 and u = 1. Consequently, the solitons of species ai experience an 
attractive force. 
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8.2 Two solitons of species OL2/CX.2 

Here we give the two-soliton solutions obtained by taking the group element E = E22 
exp [e-'"2(^i)Va2(aa2' '^i)] ^^P [6'^^'^^^''^02(^a2' -^2)]- The non-vanishing tau-functions are 

^ 1 _ g2r2(zi) _ g2r2(z2) _ ("a2 ^Q:2 + '^^2 ^^2 ) ^1 ^2 gr2(zi)+r2(z2) 



4 4 + Z2) 



2 



, "^2 '^«2 ^^2 ^^2 ^2)^ „2r2(2l)+2r2(z2) 1 1 ^ 

^ 16(^1+^2)^ ' 

_ -, _ "a2 "a2 „2r2(zi) _ ^«2 ^«2 „2r2(22) _ ("^02 ^^2 + "^0:2 ^^2 ) ^2 r2(^l)+r2(^2) 

'~ 4 ^ A ^ {zr + z^y 

I "^02 ^a2 K2 (^1 ~ ^2)"^ ^2r2(zi)+2r2(z2) ^8 19^ 

^ ^ _^ "^02 "^02 g2r2(^l) _^ ^02 ^«2 g2r2(z2) _^ ("^02 ^02 ^1 + "^02 ^^2 ^1) gr2(zi)+r2(z2) 

^ 4 4 (-Zi + -Z2)^ 

I "a2 "a2 ^a2 ^a2 (^1 ~ ^2)^ g2r2(zi)+2r2(z2) (-g ]^3^ 

16(^1 + ^2)^ ' 

^«2 _ pr2(2i) , r2(22) I Qa2 Qa2 ^a2 (^1 ~ ^2)^ 2r2(zi)+r2(z2) 

-""^e +<^a2e + 4(^1 + ^2)2 

Qga ^^2 ^^2 (^1 ~ ^2)^ -r2(zi)+2r2(22) (o 14^ 
^ 4 (., + ^2)2 ' ^^-'^^ 

_a2 _ „r2(zi) I L+ „r2(z2) I "02 ^^2 ^^2 (^1 ~ ^2) 2 r2 (zi )+r2 (22) 
-«a2e +^26 + 4{z, + Z2y 



I ^^02 ^a2 ^Q!2 



-02 ^ gA(zi) _ ^ gr2(z2) _^ '^a2 ^«2 ^«2 (^1 ^2)^ gr2(zi)+2r2(. 



^1 ^2 4^1(^1 + ^2)2 



2:2) 



_^ Qa2 "a2 ^Q2 (^1 ^2)2 g2r2(2l)+r2(22) (-g iQ^^ 



AZ2 {Zi + Z2y 

~a2 _ ^^2(21) _ ^ ^r2[z2) , "a2 ^a2 K2 (^1 ~ ^2)^ _r2(2l)+2 r2(22) 

^1 Z2 4zi(zi + ^2)2 



^ Qga "02 ^02 (^1 ^2)2 ^2r2(zi)+r2(22) (sn) 

Az2{zi + Z2y 



_ „+ ^ ^^2(21) , L+ „ „r2(z2) "02 ^«2 ^«2 ^1 (^1 ^2)2 r2(zi)+2r2(z2) 

i^,(o) - ^1 e +f^a2^2e 4(Z1 + Z2)2 ^ 

_ "a2 ^«2 ^02 (^1 ~ ^2)2 ^2 2r2(zi)+r2(z2) /'O 1 o\ 

A{z, + Z2y ' ^ ^ ^ 

_ ^a2 "a2 ^a2 (^1 ~ ^2)^ ^2 ^2 r2(2i)+r2(22) /g ^gN 

4(^1 + ^2)2 
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One can check that by taking the hmit e^^^^^' 0, such solution becomes the one-sohton 
solution (|7.5| )-(7.7). Now, considering the ratios of tau-functions appearing in the relations 

oo, such solution becomes 



( p.20| ) and ( p.27| )-( [3.32| ) one observes that in the limit e 



r2{z2) 



the one-soliton solution (|7.5| )- (|7.7| ) again, but with the replacement 



Zl - Z2 
Zl + Z2 



^.20) 



and with the ratio T2/tq changing sign. Therefore, comparing with ( |6.7| ), we observe that the 
lateral displacement and time delay for this two-soliton solution is given by ( |6.10D and (|6.11|) , 
respectively, with n = 2 and uj = 1. Consequently, the solitons of species 0^2 experience an 
attractive force. 



8.3 Two solitons of species ol^/ol^ 

Here we give the two-soliton solutions obtained by taking the group element E = = 
exp [e^^^^'^^Vag{a^,^, zi)'j exp [e^^^^^^Va^^b^,^, Z2)]. The non-vanishing tau-functions are 

_ . _ '^«3 ^"3 „2r3(^l) _ ^«3 ^03 „2r-i(z2) _ (^^3 ^^3 + '^^3 ^^3 ) r3(^l)+r3(22) 

° 4 ^ 4 ^ iz, + Z2y 

I ^^3 '^"3 ^a-i ^03 (^1-^2)"^ „2r3(2i)+2r3(22) (o 91 N 

^ iHzi + z2r ' ^ ^ 



Tl 



_^ ^a-i ^a-i g2r3(2l) ^03 ^Q3 g2r3(^2) _^ ('^03 ^^3 ^1 + "q3 ^"3 ^2) gr3(2l)+r3(^2) 

4 4 (zi + 22)2 



«a3 " ^3 ^a3 ^^3 (^1 - 2:2)^ 
16 (Zl + Z2)^ 



_^ "as "-03 "03 "03 '^^> ^2 r3{zi)+2 Tsizi) 22) 



^ ^ ]^ _^ ^Q3 ^03 g2r3(2l) _^ ^"3 ^«3 g2r3(z2) _^ ("«3 ^^3 ^1 + '^^S ^"3 ^2) ^r3{z^)+r3{z2) 

^4 4 (zi ^2)^ 

I ^^3 ^^3 ^03 ^03(^1 - ^2)^ 2r3(2i)+2r3(22) ('o 23) 

16 (zi + Z2Y ' ^ ■ ^ 



^"3 _ + pA(2l)2_, , Zoe^3(22) _ ^^3 ^^3 ^^3 ^1 ^2)^ „r3fa)+2r3(22) 

2 

8.24) 



Qga '^^3 ^^3 (^1 - ^2)^ ^2 2r3(2i)+r3(22) 



4 + Z2) 

Zl ^ Z2 " ' 42:1 {Zi + Z2Y 



-as ^ _^gr3(2i) _ ^gr3(22) _^ "as ^03 ^^3 (^1 ^2)^^ ^.r3fa)+2r3(22) 



_^ "as "«3 ^as (^1 ^2)^ g2r3(2i)+r3(22) C3 25) 
4^2(^1 + 2:2)2 



^"3 _ _„+ pr3(2i) _ r3(22) _ "«S "as ^as ^2)2 „2r3(2i)+r3(22) 



"as ^as ^as (^1 ^2)2 „rs(2i)+2r3(22) I'oofiN 
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j(^2) _ ^aa ^aa ^^3 (^1 ^2)^ „2r3(^i)+r3(z2) 



T-"3 _ _„+ 0^3(21) _ t+ r3(„^, _ - 



Qa3^a3^a3(^l r3(zi)+2r3(^2) f0 97\ 



i{z2) I "as ^03 (^1 ^2)^ ^2r3(2i)+r3(22) 



f"3 _ - pr3(2l) I L- r3(22) J 



O'asKsKA^^ ^2) r3fa)+2r3(22) (ono\ 

^"3 _ „- „r3(2i) , , - rs{z2) , "°3 "^3 ^«3 ~ ^2)^ 2r3(2l)+r3(z2) 

^iJ,(2)-«a3e +0„3e + 4(zi+^2)2 

^ Qa3 ^^3 ^«3 ~ ^2)^ ^fii(2i)+2r3(22) C3 29) 

4(2:1 + ^2)^ 

One can check that by taking the hmit e^^*^^^-* — > 0, such solution become the one-sohton 
solution ( [7.9| )-( [7?TTD . Now, considering the ratios of tau- functions appearing in the relations 



( p.20| ) and ( p.27| )- (|3.32| ) one observes that in the limit e^^'^^^) 00, such solution becomes 



the one-soliton solution (|7.9| )- ([7m| ) again, but with the replacement 



eA(.0 ^ ( ] ^m.,) (8.30) 

Zi+ Z2, 



2 



and with the ratios Ti/tq and T2/T0 changing signs. Therefore, comparing with 
observe that the lateral displacement and time delay for this two-soliton solution is given by 
( |6.10| ) and (|6.11|) , respectively, with n = 2 and a; = 1. Consequently, the solitons of species 
03 experience an attractive force. 

8.4 Two solitons of species ocxjoLi 

Here we give the two-soliton solutions obtained by taking the group element = Hyi = 
exp [e'^i^^^^VQj(a^^, 2:1)] exp [e'^2(z2)y^_^(Q±^^ ^2)] . The non-vanishing tau-functions are 



^0 



4 4 



_^ "ai "a2 "ai "^2 (^1 + ^2 ri(2i)+2 TaC^a) (g 3]^\ 

16 {zi - 2:2)2 



Ti = l I ""1 (.2 A (21) "^2 "«2 ^2r2(22) 



Q-Qi Qa2 '^Ql '^02 + ^2)^ 2ri(zi)+2r2 



16 - z^y 



e^n[zr,+2i2iz2)^ (8.32) 



_ T Qqi Qqi 2ri(2i) , Qa2 '^02 „2r2(22) 

"ai "02 '^Qi '^a2 (^1 + ^2)^ ^2ri(zi)+2r2(22) (8 33) 

16 (2:1 - 2:2)2 

f-i = a- e^^(^^) - "^i"«2 <(^i + ^2) ^r,(,,)+2r2(.2)^ (8.34) 
^ 4 (2:1 — 2:2) 



24 



:"2 _ „- „r2(22) I ^Qi ^a2 ^ai (^1 + 2 ri(zi)+r2 (za) 



L,(0) 



^^2 



'i?„(2) 



Q«i Q«2 ^1 ^2 ^ri(;.i)+r2(^2) 

- Z2 



.ai _ ^+ ri{^i) , Qa2 Qa2 (^1 + ^2) A(^i)+2r2(z2) 



a„ a 



Qa2 (^1 + ^2) ^2ri(zi)H 



_ „+ „r2(22) _ "-Qi "-a2 v^i ^ 2ri(2i)+r2(22) 

~a3 _ ^ai "'"2 „ri(^i)+r2(^2) 



- ^2 



(^1 + ^2) ri(^i)+2r2(22) 
~ai _ AUi) V ^a2 ^»2 (^1 + ^2) rifa)+2r2(^2) 

^L,{o) ~ ^ . , ^ e 



^-'-^^ Qgi Qa2 (^1 + ^2) ^2ri(^i)+r2(^2) 



-02 _ "'"2 „r2(z2) Ql "2 



^2 4 (^1 - ^2) Z2 

~a2 ^ _^gr2(22) I "ai '^02 '^ai (^1 + ^2) ^2 ri(zi)+r2(^2) 



2^2 4 (Zi - Z2) Z2 



as _ _%/^a2^ „ri(^i)+r2(22) 
Z\ — Z2 

as _ ^ai ^a2 ^2 ri(zi)+r2(z2) 
L,{2) ~ ^ 



Zl - Z2 



_ „+ ^ „A(2i) '^«2 "«2 ^1 (^1 + ^2) ri(zi)+2r2(z2) 

_ ^+ ^ „A(^i) I ^"2 Q«i '^^2 ^1 (^1 + ^2) ri(^i)+2r2(^2) 
+ 4(^1-^2) 

^"2 _ ^+ ^ „r2(^2) I ^"1 ^"1^02^2 (^1 + ^2) 2 rifa)+r2(^2) 



,r2(22) _ Q«i Q«2 ^2 (^1 + ^2) 2ri(zi)+r2(^2) 



"«-W = "-^^^ 4 (.,-.2) 

0,3 _ ^ai^a2^'^ ^ri{zi)+r2{z2) 



- Z2 



-as _ '^Q2 ^1 ri(^i)+r2(^2) 



^1 - Z2 



One can check that by taking the hmit e^^*^^^^ 0, such solution becomes the one-sohton 
solution (|7.1| )- (|7.3| ). Now, considering the ratios of tau-functions appearing in the relations 
( p.20| ) and ( p.27| )-( [3.32| ) one observes that in the limit e^'"^^^^-* oo, such solution becomes 
the one-soliton solution (|7.1| )- (|7.3| ) again, but with the replacement 

\Zl - Z2 
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and with the ratios T2/tq, t^^ /ti, t^^(-2)/t2, and t'^^q-^/tq changing signs. If we track the 
sohton a2 instead of ai, we get the same results (interchanging Fi 1^2), but the ratios of 
tau-functions changing signs are Ti/tq, t'^ jr^^ ''"r^(o)/''"0) ^^^1 t"^^(i)/ti- Therefore, comparing 
with (|6.7|) , we observe that the lateral displacement and time delay for this two-soliton solution 
is given by (|6.10|) and (|6.11| ), respectively, with n = 1 and uj = —\. Consequently, the solitons 
of species a\ and ^2 experience a repulsive force. 



8.5 Two solitons of species 0:1/0:3 

Here we give the two-soliton solutions obtained by taking the group element = 1713 
exp [e'^i^^i^Vaj (a^^, Zi)] exp \^^'^^^'^^Va^{a^^,Z'2)\. The non-vanishing tau-functions are 



^L,(0) 



4 4 

, a-^ a+ a+ {z^ - z^f 



ai "-03 "-ai "-^3 2 A (21 )+2 (22 ) 

16 (^1 + Z2Y 



7-^ = 1+ ^2 A (21) _^ ^^3 ^2 A (22) 



Q-ai '^az Q-Qi 0-03 (^1 ^2)^ g2 A (21 )+2 A (Z2 ) 
16 + ^2)=^ 



T2 = \- ""1 c^A(2i) I ^a3"^3 ^2 A (22) 



ai "-03 "-Qi "-03 V^l 2A(2i)+2 A(22) 

1 



4 

Qgi Qa3 Qa3 (^1 ~ ^2)^ 

16 (^1 + ^2)2 
_ - A(2i) I ^^3 ^^3 (^1 ~ ^2) A(2i)+2A(22) 



-as _ ^^^liWfl „A(2i)+A(22) 
^1 + ^2 



_ ^+ ^ 0^3(22) , "ai "4 (^1 ^2) ^2 2 A (21 )+ A (22) 



^ai _ ^+ pA(2i) _ ^03 '^Ql '^Q3 ^2) A(2i)+2A(22) 

^ " 4(Z1 + Z2) 



_ ^2 A(2l)+A(22) 

2^1+^2 



~a3 ^ _^gr3(22) _^ "«i "^3 "ai (^1 ^2) ^2A(2l)+A(22) 
^ ^2 42:2 (2:1 + ^2) 

~ai ^ e^l(2i) _^ '^Qi '^^3 '^Q3 (^1 ~ ^2) gA(2i)+2 A(22) 

2:1 42:1 (2:1 ^2) 

^e^l*^^^) - '^"1 '^"3 '^"3 ~ ^2) ^A(2l)+2A(22) 

^1 4^1(^1 + ^2) 



f "1 , = ^ e 



-aa ^ _^^^3_^^m_fL A(2i)+A(22) 
^2 (^1 + ^2) 
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zi + 2:2 



Wd- ^-36 4(^1 + ^2) ' ^ ' 

^"3 _ _„+ 0^3(22) , Qgi Qgi Qa3 (^^ ~ ^2) 2 ri(.i)+r3(z2) «7^ 

^ai - „+ ~, r^f^i) ^«3 ^1 (^1 ~ ^2) A(zi)+2r3(z2) ^of5o^ 

^«,(0) - «ai ^1 e 4(zi + ^2) ' ^ ^ ^ 



(zi 2:2) ^2ri(2i)+r3(22) 



^-2^^ = -2 ^A(.,)+r3(.2)^ (8_70) 



aa ^ ^^Ql^^:^lii£^^ri(2i)+r3(z2) 
^1 + Z2 

at., 



^'■^^^ -21 + 2:2 ' 

^ii,(l)-«a3e 4(Zl + Z2) ' ^ ^ 

^"3 - g- c-^«^^^^ I '^"^ ^^-^ ~ ^'^ ^2n{zi)+rz{z2) /g 73) 

One can check that by taking the hmit e^^*^^^-* — > 0, such solution become the one-sohton 
solution (|7.1|) - [7T^ . Now, considering the ratios of tau-functions appearing in the relations 
( p.20| ) and ( p.27| )- (|3.32| ) one observes that in the limit e^''^^^^-* — > oo, such solution becomes 



the one-soliton solution ( |7.1| )-( jr3| ) again, but with the replacement 



e^U^iJ ^ 11 — 11 e'll-i) (8.74) 

and with the ratios Ti/tq, T2/ro, r^'^ /ti, t'^^(^2)/'^^^ changing signs. If we track the 

soliton ^3 instead of ai, we get the same results (interchanging Fi ^ F^), but the ratios of tau- 
functions changing signs are Ti/tq, t'^ /tq, and t''^*^2)I'^'2- Therefore, comparing with 
( |6.7|) , we observe that the lateral displacement and time delay for this two-soliton solution is 
given by ( |6.1CI| ) and ( |6.11| ), respectively, with n = 1 and u = 1. Consequently, the solitons of 



species cti and cts experience an attractive force. 



8.6 Two solitons of species 0:2/0:3 



Here we give the two-soliton solutions obtained by taking the group element S = E23 
exp [e''"2^^^-*V'Q,2(a^25 ^1)] ^xp [e^^^^^^Vas{a'^^, Z2)]- The non-vanishing tau-functions are 



_ _ Q-aa ^a2 J2 Tiiz^) _ ^03 "«3 „2r3(22) 

^° 4 ^ 4 ^ 



Qga ^as ^^2 ^as (^1 ^2)^ ^2 r2(2i)+2 r3(22) (8 75) 

16 {Zi + Z2)'^ 
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4 4 



Qa2 "aa '^^2 '^03 ^2)^ g2r2(zi)+2r3(22) 



16 (zi + ^2)2 
4 " '4 



-j^ _^ Qa2 Qa2 g2r2(2i) _^ "aa "aa g2r3(z2) 



^a2 "^aa '^02 ^^3 (^1 ^2)^ 2 r2(2i)+2 r3(z2) 

16 {zi + z^Y 



""2 ^1 „r2(2l)+r3(z2) 



-Zl + Z2 



^012 ^aa ^a3 (^1 ~ ^2) r2(zi)+2r3(22) 



4 (2:1 + ^2) 

„+ y ^r3(z2) I "q2 "q2 '^aa ~ Z2 2r2(zi)+r3(z2) 

' 4{z^ + Z2) 

Qq2 ^X, ^2 r2(zi)+r3(22) 

zi + 2:2 

p^2(^l) _ "aa "a2 (^1 ~ ^2) r2(^i)+2r3(^2) 

4(Z1 + Z2) 

"«a gr3(^2) _^ '^«2 "^aa '^«2 (^1 ~ ^2) ^2 r2(^i)+r3(z2) 

^2 4 2;2 (2;i + 2:2) 

^2 (^1 + ^2) 



^"3 "q2 --^ „r2{zi)+r3{z2) 
^ 1 



"q2 „r2(2l)+r3fe) 

*5 ) 



-Zi + Z2 

gr2(zi) _^ Q«2 "«a "aa (^1 ~ ^2) gr2(2i)+2 r3(z2) 
Zi 4:Zi{zi+Z2) 

_^gr2(zi) _ ^a2 ^O'i °«a ~ ^2) ^r2(2l)+2r3(z2) 

2^1 4^1(^1 + ^2) 

r^3(z2) I '^«2 '^«2 "^aa (^^ ~ ^2) 2 r2(^i)+r3(^2) 
^ 4(^1 + ^2) 

^^3(22) _ "a2 ^02 "^^3 (^1 ~ ^2) 2r2(2i)+r3(22) 
4(^1 + ^2) 

Qga Qa3 ^1 ^2 r2(zi)+r3(z2) 

Zl + Z2 
— + 2 
Qga "aa ^2 _r2(zi)+r3(z2) 

Zl + Z2 

„+ ^ ^r2(zi) "aa "02 "^aa ^1 (^1 ~ ^2) r2(zi)+2 r3(22) 

' HZ, + Z2) 



,r2(zi) _^ ^aa '^a2 '^^a (^1 ^2) ^r2(0i)+2 TaCzj) 

4 (^1 + ^2) 
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fla, O-lo \Zl - Z2) 



^"3 _ „- p-ns(22) I "-«2 "-^3 "^02 
- "a3 ^ 



,2r2(^i)+r3(22) 



4(;2i + Z2) 

r3(^2) _ "a2 "a3 "^2 (^1 ~ 2r2(^i)+r3(^2) 



4 Ui + Z2\ 



(8.94) 
(8.95) 



One can check that by taking the hmit e^^*-^^-* — > 0, such solution become the one-sohton 
solution (|7.5| )- ([7T7| ). Now, considering the ratios of tau-functions appearing in the relations 
( p.20| ) and ( p.27| )- (|3.32| ) one observes that in the limit e^'"^*^^^^ 00, such solution becomes 



the one-soliton solution (|7.5|)-([7I7|) again, but with the replacement 



.^2(^1) 



Zl - Z2 
Zl + Z2 



^.96) 



and with the ratios Ti/tq, T2/T0, t2'^/t2, and 't^^o^/tq- If we track the soliton 03 

instead of a2, we get the same results (interchanging /2 ^ -^3), but the ratios of tau-functions 
changing signs are T2/TQ, t^^/tq, and 'r^%)/T2- Therefore, comparing with 



we 



observe that the lateral displacement and time delay for this two-soliton solution is given by 
( |6.10| ) and (|6.11|) , respectively, with n = 1 and u = 1. Consequently, the solitons of species 
a2 and 0:3 experience an attractive force. 



8.7 Summary of time delay results 

The calculations of the time delays of the two-soliton solutions have shown that the solitons 
of the same species experience an attractive force. The pair of solitons of species ai and 0:3, 
as well as 02 and also experience attractive forces. However, the pair of solitons ai and 
a2 suffer a repulsive force. In addition, the force between sohtons of the same species is twice 
as strong than that between solitons of different species. Consequently, the sign u and the 
integer n appearing in the time delay ( |6.11 ) is determined by the scalar product of the roots 
by the formula {ai\aj) = cun, by normalizing the roots as af = 2. Therefore, the time delay 
suffered by a sohton of species a,, with rapidity 9i when scattering with a soliton of species 

> 6'2, is given by the formula 



aj, and rapidity 62, with 9 



{ai\aj) 
rrii cosh 61 



In 



tanh 



61 — 62 



(8.97) 
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A The affine Kac-Moody algebra s 13(C) 

In this appendix we give the necessary information about the affine Kac-Moody algebra 



s[3(C). More details can be found in the book by Kac which we follow in our presentation 
of basic facts and definitions, or in the review paper by Goddard and Olive intended for 
physicists. 

Recall that the Lie algebra sl3(C) is the Lie algebra formed by all 3 x 3 complex matrices 
with zero trace. We use the standard choice for a Cartan subalgebra () and for a basis of 
the corresponding root system A. Denote the three positive roots by a,, i = 1,2, 3, with 
a = 1, 2, being the simple roots and 03 = ai + 02- The Cartan generators are 





















-1 










1 




V 





0^ 




^ 








ifi = -1 , H2= \ Q I , (A.l) 
\0 00/ \00-l/ 

and the basis vectors of the root subspaces corresponding to the positive roots are chosen as 

/010\ /000\ /00l\ 

E+^, = 0, = 001, E+„3 = 0. (A.2) 






For negative roots we have 

= {E+^f. (A.3) 
We use the invariant bilinear form on s[3(C) defined by 

{x\y)=ti{xy), x,yesh{C). (A.4) 

The restriction of this form to the Cartan subalgebra 1) is nondegenerate, therefore it induces 
a nondegenerate bilinear form on f)* which we denote by (■ | ■) as well. Note that with this 
definition one has 

(ai|Q;i)=2, (a2|a2) = 2, (ai|a2) = -1- (A.5) 

The affine Kac-Moody algebra shiC) is constructed in the following way. Consider the 
loop algebra 

£(s[3(C)) = C[C,C-']®5[3(C), (A.6) 

where C[C, C~^] is the algebra of Laurent polynomials in C^. An element of i2(5[3(C)) is a finite 
linear combination of the elements of the form x, where m G Z and x G sl^{C). The 
structure of a Lie algebra in C{s[^{G)) is introduced by the relation 

[r®x,C®y] = C^''®[x,y]. (A.7) 

We identify the Lie algebra s[^{'C) with the subalgebra of C{s{^{'C)) formed by the elements 
of the form 1 ® x, a; G s[3(C). This allows us to write CJ^x instead of ® x. Similar 
identifications are used below. 

Define a C- valued 2-cocycle on C{s[^{£-)) by 

V'lCx, Cy)=m(x|i/)W«,o (A.8) 



30 



and denote by £(s[3(C)) the extension of C^sl^^C)) by a one-dimensional center associated 
to the cocycle ip. Denote by C the corresponding center element, then the commutation 
relations in the Lie algebra £(s[3(C)) = C^sl^^C)) © CC are given by 

[Cx, Cy] = r'-'^ix, y]+mC6m-,n,o- (A.9) 

Now denote by sl3(C) the Lie algebra which is obtained by adjoining to £(s[3(C)) a derivation 

D = ({d/d(). The commutation relations for the Lie algebra sl^lC) are defined by relations 
( |A.9| ) and by the equalities 

[D,Cx]=mCx, [D,C] = 0. (A.IO) 

The elements 

H^ = rH,, H^ = rH2, E^ = rEo., (A.H) 

where m G Z and a E A, together with the elements C and D form a basis for s[3(C). These 
elements satisfy the commutation relations 

[i/r, H^]=m{aa\ai,)C6m+n,0: (A.12) 

[^r,^±„J = ±(«a|a.)i?rar, (A.13) 

[EZ. E\J = Hr- + mC5^^^,o. (A.14) 
[E^, E\,^ ] = //r™ + i/r"-" + mC5^+.,o, (A.15) 

[^rp^:J = ^^3+^ (A.16) 

[E-,E!!,J = -E-+", (A.17) 
[EZ.E\^] = E^^\ (A.18) 

where a, 6 = 1,2, i = 1,2,3 and m,n G Z. The remaining non-vanishing commutation 
relations are obtained by using relation ( |A.3| ) which implies that 

[^r,^;] = -[^-.,^-/3]- (A.19) 

We use in the paper the principal Z-gradation of s[3(C) which can be defined by the 
grading operator 

Qppai = H^ + H2 + ^D (A.20) 
in the following way. Consider the following subspaces of s[3(C): 

0m = {x e s[3(C) I [Qppai, x] = mx}, m G Z. (A. 21) 

It is easy to get convinced that 

55(C) = 0S™, (A.22) 

and, moreover, 

[flm, 0n] C 0m+„. (A. 23) 

Hence, we really have a Z-gradation of 5l^{<C). The subspace go is a actually a subalgebra of 
s[3(C) given by 

flo = C/fi©C/f2©CC©CD = C/7i©Cif2©CC©Cgppai, (A.24) 
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and for other subspaces Qm we have 

Qsm = CH^(BC H^, m^O, (A.25) 
S3„.+i = CE:^^®CE^^®C , (A.26) 
gsrr.+2 = C © C ®^Ks- (A.27) 

Among the representations of g there are three which play an important role. They are 
the fundamental highest weight representations. Each one has a highest weight state \Xj), 
J = 0, 1, 2, satisfying 

Ha\Xo)=0, Ha\Xb) = 6a,b\\), C\X,) = \X,) (A.28) 

for a, 6 = 1, 2 and j = 0, 1, 2. Such states are annihilated by all positive grade subspaces 

gm\Xj)=0, m> 0, (A.29) 

and all the states of the representation spaces are obtained by acting on the corresponding 
highest weight state | Aj ), with negative grade generators. The representation spaces can be 
supplied with a scalar product with respect to which 

(i/r)^ = H-^. = EZ^a, (A.30) 

Ct = C, = D. (A.31) 

It follows from (|A.24|) - (|A.27|) that in this case 

(gm)^ = S-n^, (A.32) 

and one has 

( A, I = 0, m > 0. (A.33) 

Recall that go is a subalgebra of sl^iC). It is convenient also to consider two additional 
subalgebras 

0<o = 0m, 0>o = 0-m- (A.34) 

m>0 m>0 

These subalgebras and the corresponding Lie groups play important role in the dressing 
transformation method. 
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B Mathematica program implementing Hirota method 



In this appendix we give the source of a Mathematica program which allow to obtain and 
test soliton solutions of equations (p.l|)-( p3D using the Hirota method. 

First of all note that to use it one has fix the number of solitons and the soliton species. 
This information is saved in the variable n and in the list s respectively. The Hirota ansatz for 
tau- functions is given by the function t. The corresponding coefficents are collected into the 
arrays b. The function eqs extracts from the equations for the tau-functions the terms with 
the fixed total power of The function calcbs calculates the coefficients of the Hirota 
ansatz iteratively. The function testsol test that the equations for the tau-functions are 
satified by the solution found. The concrete code follows. 

g[i_, z_, xp_, xm_] := m[i] (z xp - xm/z)/2; 
m[3] = m[l] + m[2] ; 

c[l__][i_] := 1 /; i <= 3 && Plus [1] == 0; 
c[l__] := /; i > 3 && Plus [1] == 0; 

Ind[p_] := Select [Allind, Plus 00 #1 == p &] ; 
b[l__] := Array [c[l] , 21] ; 
t [k_ , xp_ , xm_] : = 

Sum [Product [(e[i] Exp [g [s [ [i] ] , z [i] , xp, xm] ]) "Allind [ [j] ][ [i] ] , 

{i, n}](c 00 Allind [[j]]) [k] , { j , Length [Allind] >] ; 
eqs[p_, 1_] : = 

Block [{temp = Ind [p] [ [1] ] } , 
Table [D 00 

Prepend [Table [{e [i] , temp[[i]]}, {i, n}] , 

taueq[k, xp, xm]] /. {e [i_] -> 0}, {k, 21}]]; 

calcbs := 

Block [{temp, outf}. 
Off [Solve : : svars] ; 
Allind = Flatten [ 

Table 00 Join [{Table [i [k] , {k, n}]}. Table [{i [k] , 0, 2>, {k, n}]], 
n - 1] ; 

Do[ 

temp = b 00 Ind[p] [[k]] ; 
Evaluate [b 00 Ind[p][[k]]] = 

Factor [temp /. Solve [eqs [p, k] ==0, temp][[l]]], {p, 2 n}, 
{k. Length [Ind [p] ]}] ; 
outf = ToString[n] <> "sol" <> ToString [SequenceForm 00 s] <> ".m"; 
If [Length [FileNames [outf]] != 0, DeleteFile [outf ] , Null]; 
Save [outf , c] ; 

Print ["The results of calculations were saved in file " <> outf]; 
On [Solve : : svars] ] ; 
testsol := 

If [Table [Simplify [taueq[k, xp, xm] /. {xm -> 0, xp -> 0}], {k, 21}] 
== Table [0, {21}] , 
Print ["Equations are satisfied"]. 
Print ["Equations are not satisfied"]]; 

Below we give the script of a sample Mathematica session. It is supposed that the above 
code is in the file gensol .m and the equations for the tau-functions are in the file taueqs .m. 
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In[l] := « taueqs.m 
In [2] := « gensol.m 
In[3] := n = 2; 
In[4] := s = {1, 3}; 

In [5] := Calebs 

The results of calculations were saved in file 2soll3.m 

In [6] := testsol 
Equations are satisfied 

In conclusion note that the program above works in principle for any number of solitons, but 
already for three-soliton solutions the procedure of testing the solution is very time consuming. 
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